UTAS-PHYS-2007-22 
December 2007 



The molecular asymmetric rigid rotor Hamiltonian as an exactly 

solvable model 



P. D. JarvitS and L. A. Yates 

School of Mathematics and Physics 
University of Tasmania, Private Bag 37 
7001 Hobart, Tasmania, Australia 
E-mail: Peter . JarvisSutas . edu . au ; Luke . YatesSutas . edu . au 

Abstract 

Representations of the rotation group may be formulated in second-quantised language via 
Schwinger's transcription of angular momentum states onto states of an effective two-dimensional 
oscillator. In the case of the molecular asymmetric rigid rotor, by projecting onto the state space 
of rigid body rotations, the standard Ray Hamiltonian H(l, K, —1) (with asymmetry parameter 
1 > k > — 1), becomes a quadratic polynomial in the generators of the associated dynamical 
su(l, 1) algebra. We point out that H(1,k, — 1) is in fact quadratic in the Gaudin operators 
arising from the quasiclassical limit of an associated su q (l, 1) Yang-Baxter algebra. The general 
asymmetric rigid rotor Hamiltonian is thus an exactly solvable model. This fact has important 
implications for the structure of the spectrum, as well as for the eigenstates and correlation 
functions of the model. 
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1 Introduction 



The analysis of rotational dynamics is central to molecular spectroscopy and has been much studied 
in the literature. The paradigmatic model is of course the rigid rotor, whose symmetric (prolate- 
oblate) limit is trivially handled in the quantum case using the angular momentum algebra. For 
the asymmetric case one has an additional parameter k, and although, in principle, the model is 
easily solved for given angular momentum by a matrix diagonalisation [H[2]) it is surprisingly rich. 
Classical approaches to the full, non-rigid case achieve a separation of vibrational and rotational 
degrees of freedom by imposing subtle geometrical constraints on the relative coordinates, which 
often generate expansions of considerable combinatorial complexity when accurate descriptions of 
couplings are developed O IU [5]. Such is the fundamental importance of the general problem 
that new contributions continue to appear, as for example the recent work [6] concerning large N 
expansions, and [7] on a semiclassical approach to expansions of the energy levels and wavefunctions 
near the symmetric limits k = ±1. 

In contrast to ab initio approaches are algebraic developments which adopt as a starting 
point, an appropriate Lie algebra of operators with enough structure to capture important features 
of the molecular Hamiltonian (such as interactions which dominate near level degeneracies), but 
where results from representation theory can be a valuable organising framework. An early study of 
potential dynamical symmetry groups for the rotor problem is jSj; unitary dynamical symmetries 
have been extensively used in [9]. In the end such treatments, as well as the above-mentioned 
formal expansions from classical dynamics, tend to become semi-empirical studies of appropriate 
parametrisations in an attempt to reproduce specific aspects of observed data. 

A problem equivalent to the rigid rotor model was studied some time ago by Patera and 
Winternitz [10] in the context of rotation group representation basis labelling. In that work, the 
system was shown to be separable in elliptic coordinates on the sphere, and the eigenfunctions and 
eigenvalues turned out to be related to certain Lame polynomials, with the asymmetry k being 
essentially the modulus for the Jacobian elliptic functions involved in the coordinate system. Thus, 
to the extent that the relevant Lame polynomials are available and tabulated, the problem can be 
regarded as completely solved in principle by this analytical method, in a similar way that it can 
be said to be solved on the algebraic side, by a straightforward matrix diagonalisation. 

The present work is intermediate between the classical dynamics-inspired works, and the use 
of algebraic models as dynamical symmetries, while complementing the analytical approach. It 
extends the well-known, textbook fact that the symmetric (oblate-prolate) limit of the rigid rotor 
is diagonalisable in standard angular momentum bases, to the class of molecular asymmetric rotor 
Hamiltonians, which are shown to be exactly solvable in the sense of Bethe Ansatz techniques (see 
below). 

In ^2] below we review some transcriptions of the rotor problem using well-known algebraic 
tools, chiefly the model of Schwinger [TT| which gives a realisation of the su(2) angular momentum 
algebra in the second-quantised language of modes of a two-dimensional oscillator. The relation to 
the algebra of the non-compact symmetry group su(l, 1) is also noted [12]. In Sj3] the symmetric and 
asymmetric rigid rotor Hamiltonians (in the Ray [I] parametrisation H(l,n,—1) ) are compared 
with the class of so-called Gaudin operators [13J. These are terms in the expansion of transfer 
operators coming from a specific realisation of a certain Yang-Baxter algebra in the so-called quasi- 
classical limit. To the extent that they are coefficients in u expansions of formal operators t(u) 
which mutually commute, [t(u),t(v)} = for all u,v, the Gaudin operators, or any polynomials 
thereof, are simultaneously diagonalisable via Bethe Ansatz techniques in the tradition of exactly- 
solvable models in statistical and condensed matter physics [14 } 115 1 fT6] . The symmetric rotor (as 
might be expected) can trivially be identified with such transfer matrix coefficients for the two- 
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dimensional Yang-Baxter algebra with rational i?-matrix (the su(l, 1) case), whereas the extension 
to the trigonometric su q (l, 1) case yields, in the quasi-classical Gaudin limit, the asymmetric rigid 
rotor Hamiltonian. The deformation q = e 1 " 1 provides the additional asymmetry parameter, the 
identification being essentially cos 7 = k (after rescaling 7 — > 777, u — > u/r/ and taking the limit 

The conclusion of our investigation is therefore that the molecular asymmetric rigid rotor does 
indeed belong to the exactly solvable class. This fact has important implications for the structure 
of the spectrum, as well as for the eigenstates and correlation functions of the model. Our work also 
establishes hitherto unnoticed connections between the analytical approach via special functions, 
and integrable models of this type. Further discussion and conclusions are given in £JH below. 



2 Oscillator formalism 

The energy operator for the asymmetric rigid rotor when written in body-fixed coordinates adapted 
to the principal moments of inertia of the rigid body takes the form 

H = + — + 

2I a 21b 2I C 

where J a ,Jb,J c are the generators of rotations about the body-fixed axes, and I a ,Ib,I c are the 
principal moments of inertia. For future reference we define a = (2J tt ) , b = (2J&) , c = (2I C )~ 1 , 
so that the Hamiltonian becomes 

H(a, b, c) = aJ 2 a + 6J 6 2 + cJ c 2 . (1) 

The standard way of writing the rigid rotor Hamiltonian is via the so-called Ray [Tj Hamiltonian 
H(1,k,—1) defined by deriving the simple relation 



H(a, b, c) = i(a - c)H(l, k, -1) + ±(a + c)H(l, 1, 1), (2) 
2b- a- c 

K = , 1 > K > — 1 



where it is assumed that a > b > c, and a fortiori, 1 > k > —1, and H(l, 1, 1) plays the role of 
the Casimir operator C{su{2)). In this case the oblate and prolate cases are k = +1 and k = — 1, 
respectively. 

States of the system are spanned by the set \J, M.; j, m), where J {J + 1) is the eigenvalue 
of the square of the total (centre-of-mass) angular momentum operator, with J the total angular 
momentum quantum number, and Ai its third component; we have as usual M = —J\ — J + 
1, • • • , J — 1, J. j = J is similarly the angular momentum quantum number for body-fixed 
rotations. Thus each eigenstate of H(a, b, c) is 2J + 1-fold degenerate, corresponding to the values 
of «M, and the Hamiltonian must be diagonalised, for fixed J = j, identically for each A4, only in 
the (2j + 1) x (2j + 1) subspace of states of differing m, m = —j, — j + 1, • • • — We henceforth 
drop explicit reference to the J ', M labels and refer to the basis simply as | j, m) - in other words, 
we are exploiting conservation of angular momentum, and referring the solution of the complete 
rigid rotor problem on the space ©j =J=0 Hj®Hj, to the diagonalisation of an effective rigid rotor 
Hamiltonian on the model space (B^qHj . The energy eigenstates can be denoted \ \j, r)), where r 
is a pseudo-magnetic quantum number labelling states of increasing energy and running over 2j + 1 
values —j, — j + 1, • • • , j (see [HE]). In this notation it is known PQ that the corresponding energy 
eigenvalues E K (j,r) satisfy the reflection property, E K (j,T) = E^ K (j, — r), so that the rigid rotor 



2 



problem essentially amounts to diagonalising the one-parameter family H(l, k, — 1) in the restricted 
interval < K < 1. 

In the formulation of Schwinger [11], the su(2) algebra of angular momentum can be expressed 
in terms of the raising and lowering modes of a set of two independent oscillators (which can be 
thought of as associated to an auxiliary two-dimensional space). Calling these mode operators a, , 
b, , the angular momentum generators are 

J x = - iiffib - b^a) , J y = \ (a>b + V a) , J z = \ (ate - b^b) (3) 

which can be easily checked to fulfil the required su(2) Lie algebra relations, 

[J x , J y ] = iJ z , and cyclically, 
provided the oscillator raising and lowering operators fulfil the usual relations 

\a,a'] =1, [a, a] = [a^a^] = 0, 

[b,V]=i, [lb] = [V,V} = o, 

with of course 

[a,b]=[a,V] = [a\b] = [a\V]=0. (4) 

The advantage of the Schwinger transcription is that angular momentum eigenstates \j, m) 
as above can explicitly be written down in terms of Fock mode states created by monomials in the 
a) , W operators acting on the vacuum state |0, 0} defined by 

a|0,0) = 6|0,0) = 0; 
(a) ) Na &) Nb 

\N A ,N B ) =i-2_lJ_|0,0), N A ,N B = 0,1,2,- ■■ , 

with N A \N A ,N B ) =N A \N A ,N B ), N b \N a ,N b ) = N b \N a ,N b ), where 

N A :=a ] a, N B :=Vb. (5) 

Moreover, it is easily seen that the angular momentum quantum number itself is intimately related 
to the total number operator for oscillator modes 

N=N A + N B . (6) 

Explicitly, an easy calculation gives 

C(su(2)) = Jl + J 2 y + Jl = \N{\N + 1), (7) 

so that the angular momentum operator is j = |iV with eigenvalue j = = ^(^,4 + Nb) on the 
above states. Thus, in terms of the \ j, m) states we have in fact 

\N A , N B ) =\j = \{N A + N B ), m = \(N A - N B ) ). (8) 

Any operator acting on states of the rotor system can now be expressed in terms of monomials in 
the creation and and annihilation operators in the equivalent number basis, the asymmetric rotor 
Hamiltonian itself being a case in point. 

Before proceeding, it is useful to make a different operator transcription, this time from 
the Schwinger angular momentum generators to the generators of a closely related algebra, also 
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realised in terms of the oscillator modes, namely the Lie algebra su(l,l) [12] . The latter algebra 
has generators K + , , Kq and nonzero commutation relations 

[K+, KJ[ = - 2K , [K , K±] = ±K±, 

which are fulfilled by the following operators acting on the A, B Fock spaces separately: 

K A = iStat, K A = |ao, K A = I^a + ±, 

K B = ±W, K B = ±bb, K B = ^b + \, (9) 

as can easily be seen from the fundamental commutation relations of the oscillator modes. The 
Casimir operator 

2C(«u(l, 1)) = K + K_ + K^K + - 2Kl 

for the representations of su(l, 1) on each oscillator space has fixed eigenvalue C A = — ^ = C B 
(corresponding to the common eigenvalue k(k + 1) for the direct sum of two irreducible representa- 
tions with spin quantum numbers k = — j, — § in each case). By contrast, the total su(l, 1) algebra 
with generators K± = K A + K±, Kq = Kq + Kq has Casimir 

C A+B (su(l, 1)) = l(N A - N B f - \ (10) 

(corresponding to spin quantum number k A+B = — ^iV/i — A^^l — ^)- 

From the above definitions we finally have general second-quantised forms for the asymmetric 
rigid rotor Hamiltonian, once an identification between J x ,Jy,Jz an d J a , J&, J c is made, which 
is to say an identification between the x,y,z directions and the (oriented) principal body axes 
corresponding to a,b,c. The canonical choice, in view of the distinguished role of the b axis (in the 
Ray parametrisation) vis a vis the z axis (with respect to the standard choice of deformation of 
the su{2) Lie algebra, as will be seen below), turns out to be J x = J a , J y = —J c , and J z = J&, 
giving straightforwardly 

H(a, b, c) = \{a - c) ( J a 2 - J c 2 ) + \ (a + c) ( J a 2 + J c 2 ) + b J 2 

= (a - c)(K A K B + K A K B ) + (a + c){2K*K* - J) + b{K* - K B f 

= (a - c)((K A K B + K A K B ) - 2kK a K b ) + b(K A + K B f (11) 

(up to a constant), where the last rearrangement also presages the discussion to follow. 

3 Exactly solvable rotor Hamiltonians from su(l, 1) and su q (l,l) 
Gaudin operators. 

In the theory of integrable quantum Hamiltonians |14t \TE[ [T6] , a family of simultaneously diago- 
nalisable operators may be generated from a set of algebraic equations called the RTT relations. 
Exactly solvable Hamiltonians may be constructed as polynomials in these operators. The struc- 
ture of the relations is determined by a given .R-matrix, a (numerical) solution to the Yang-Baxter 
equation, and the RTT relations themselves generate the Yang-Baxter algebra. Here we consider 
for the two dimensional case the so called rational solution, its generalisation to the trigonometric 
solution, which for appropriate real forms generate the algebras Y(su(l, 1)) and Y(su q (l, 1)) re- 
spectively, and finally the so-called quasiclassical limit of the latter. Details of the constructions 
are given in the reviews cited; see for example also |17j. 
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3.1 Rational i?-matrix: su(l, 1) transfer operator and the symmetric rotor. 



The rational solution of the Yang-Baxter equation in the two-dimensional case is the 4x4 numerical 
R matrix, 



R(u) = - 
u 



f u + T] \ 

u rj 

T] u 

\ u + rj j 



(12) 

which generates the Yang-Baxter algebra via the RTT relation using the 2 x 2 L operator 





u + r) J 


7]J^ 






U - 7]J 




U + T]Kq 




:( 


, -VK+ 


U — T]K{ 



L K (u) =- ( " ' ^ ) , (13) 

u\ -T]K + u - 7]K J 

where in the compact case J±, Jq are generators of the su{2) Lie algebra, and in the noncompact 
case the K$,K± generate su(l,l). In either case the monodromy operators T(u) are obtained as 
matrix products; for the present case we have simply 

T(u) = ( ^ e \ ) • L A (u - e A ) ■ L B (u - e B ) 

for a tensor product of just two copies A, B of the generators, where 5,ea,£b are parameters; the 
above expressions are explicitly functions of u/rj, in order to facilitate expansion of the operators 
about u — > oo via the limit r/ — ► of the parameter r\. 

From (I12p . (|13p we have for the transfer operator, the trace t(u) = tr(T(u)), 
t(u) = [2 + r 1 2 5 2 ]+r 1 2 5[-^-+ K * 1 



We define 



obtaining 



-(u-e A ) {u-e B ) 



T A = lmi t(u), 

(u — £a)(u — £b) , x -n 
r B = lim i ^ —t{u), (14) 



r A = 5(e A - e B )K$ + 2K*K B - (K^K B + tfftff), 

t b = - S(e A - e B )K B + 2K*K B - {K*K B + K*K B ). (15) 



In comparison with the asymmetric rotor Hamiltonian, setting Ae = {ea — £b) as the (arbitrary) 
spectral parameter shift, it is clear that the combinations \{ta + tb) and (ja — t b)/(2SA£), will 
produce (after setting the parameter 5 to 0) the independent bracketed terms in (jlip . with however 
the restriction for the first, Casimir-type term, that n = +1; that is, the oblate a = b > c case. (The 
prolate case follows with use of the general reflection property mentioned above) . The conclusion of 
this reformulation, as might be expected for this rather trivial, non-deformed case, is that the rigid 
symmetric (prolate/oblate) rotor is recovered as an instance of an exactly-solvable model associated 
with the su(l, 1) Yang-Baxter algebra. 
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3.2 Trigonometric i?-matrix: su q (l, 1) transfer operator, quasiclassical limit, and 
the asymmetric rotor. 

To go further using the same algebraic template, an additional parameter is needed. This is 
provided by the trigonometric generalisation of the rational solution of the Yang-Baxter equation. 
The trigonometric i?-matrix which generalises (|12p . is the numerical 4x4 matrix [15] 



R(u; 7 ) 



1 



sin(7u) 



/ sin(7(n + r ? )) \ 

sin(7u) sin(7?7) 

sin(7ry) sin(7n) 

\ sin(7(n + ??)) / 



which now is a function of the (scaled) quantities u/r] and the deformation q = e* 7r? . The corre- 
sponding L-operator is 



L J (u; 7 ) 



sin(7(-u + t]Jq)) sin(7ry)J_ 



sin(7«) I sin(7r/)J + sin(7(u — rjJo)) 



where Jo, J± satisfy the su q {2) algebra 



[J ,J±] = ±J±, [J+,J-] 



sin(27J ) 
sin 7 



or alternatively 



sin(7ti) I 

where Kq,K± satisfy the su q (l, 1) algebra 

[K ,K±] = ±K±, [K+,K. 



sin(7(it + rjKo)) sin^^K- 
—sin( r yr])K + sin(7(u — T]Kq)) 



sin 7 



We note the expansion, 



L(u;i) 



1 
1 



+ rj 



7 / cos(7«)l£o 



(16) 



sin(7«) 



-if 4 



K_ 
cos(7«)ii"o 



Wl 2 {K f 



1 
1 



0{if 



so that in addition to the standard asymptotic limits of L (and similarly R) as unit matrices for 
u —* 00 (or Tj — > 0), we see that for 7 —* (q —* 1) the trigonometric solutions revert to the ordinary 
rational case: 



hm L(u; 7) = - ' 

7-»0 It V — r l^+ u — V^0 



(17) 



Now consider the following monodromy operator T(u) defined as a 2 x 2 matrix product over the 
enveloping algebra of the K A and K B generators, 



T(u) 



e Sr > 
e - *» 



L A (u-£ A )-L B {u-e B ) 
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(suppressing the 7 label in the notation). For present purposes it is sufficient to record merely the 
transfer matrix t(u) := tr(T(u)) in the scaling limit 77 — > 0, 

t(u) = [2 + V 2 5 2 ] + 2r ? V[cot( 7 (u - e^))^ + cot( 7 (u - e B ))K B ] 

+ r/V [2cot( 7 (n - e A )) cot( 7 (n - s b ))k£k^ - csc( 7 (-u - e A )) csc( 7 (u - £b))(^^ B + )] 
-f[[Ktf + {K§)*\+Otf)i 

note that we recover the standard su(l, 1) Lie algebra from (|16p as g — > 1, or just 77 — >• 0. We define 

sin( 7 (u-e yl ))sin( 7 (n-eB)) 
ta = hm ^ tlu), 

y sm(j(u - e A )) sin( 7 (u - eg)) 
t b = hm 5 tlit), (18) 

obtaining 

t a =25 1 sm( 1 (e A -e B ))K£ + 2 7 2 cos^a-^))^^ " l\K^K B +K^K^), 

t b = -25 7 sin( 7 (eA-^))^o B + 2 7 2 cos( 7 (eB-e A ))^o A ^o iJ " ~f 2 (K^K B + K A K B ). (19) 

Finally we can compare appropriate combinations of ta, tb with the desired asymmetric rigid rotor 
in second-quantised form. Again \{ta + tb) / 7 2 and (ta — rs)/(2 7 <5Ae) provide the correct material 
(after setting the parameter 5 to 0) to produce the two bracketed terms in (fTTj) ; in contrast to the 
situation for (|15p however, there is now no longer a restriction on the relative weights within the 
first, Casimir-type term, and we identify 

k := cos( 7 Ae). (20) 

Thus, the molecular asymmetric rigid rotor Hamiltonian can be recovered as a polynomial in the 
independent Gaudin operators t a ,tb- The deformation parameter 7, in the quasiclassical scaling 
limit, is thereby given via (I20p in terms of the Ray asymmetry parameter k, in units defined by 
the spectral parameter shift Ae which plays the role of an (arbitrary) auxiliary parameter in the 
reformulation. 



4 Discussion 



In this note we have established that the classic molecular asymmetric rigid rotor Hamiltonian be- 
longs to the exactly-solvable class of models, in the tradition of statistical and many body condensed 
matter physics. As pointed out in the introduction, the problem is also 'solvable' in principle, in 
the usual sense, either algebraically via matrix diagonalisation, or analytically via special functions 
as solution families of differential equations deriving from Schrodinger's equation. These standard 
methods of course increase in complexity with increasing angular momentum, and do not have 
much to say about the general structure of the eigenfunctions and eigenvalues. 

By contrast, our reformulation of the model as an exactly solvable system has important 
implications for the general structure of the solutions. For example, the algebraic Bethe Ansatz 
implies that the eigenfunctions (for arbitrary angular momentum) can be given a direct product 
form, with parameters satisfying a set of algebraic equations (in the quasiclassical Gaudin limit); 
the energy eigenvalues are in turn able to be written down explicitly as polynomials in these 
Ansatz parameters. Moreover, the Bethe Ansatz solutions also open the way to exact evaluation 
of correlation functions, and hence to the physical properties of the system. At the same time, our 
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work implies that there are close and hitherto unnoticed links between the manipulations involved 
in the special function methods applicable in this case (such as continued fraction expansions for 
the Lame functions, as discussed in [TO]), and exactly solvable models. 

At this point our analysis is on a par with analogous studies pointing out the relevance of 
exactly-solvable models in a range of problems in many-body and condensed-matter contexts (see 
for example the review [17] and the paper [18] on Gaudin algebras). One point to note is that, 
whereas the iV-body systems often involve spin chains with a number of systems with finite state 
spaces, (and for example, su(2) symmetry algebras), in our second-quantised reformulation, the 
'chain' consists of just two systems A and B, but with each having an infinite number of states. 
Moreover, the relevant symmetry algebra is the non-compact su(l, 1) Lie algebra, in contrast to 
what might have been expected, given the obvious analogy between the asymmetric rigid rotor 
and anisotropic 'XXZ' or 'XYZ' type spin chains where deformations of su(2) are involved. For 
further details of the notation and main results needed from exactly solvable models in many body 
quantum theory, we refer the reader the literature. 

One final comment should be made regarding the status of the work in the context of group 
representation theory, which was the focus of the paper [TO]. In that paper the central question was 
the use of a non-subgroup labelling scheme for representations of the rotation group SO (3) = SU (2). 
Our transcription shows that the choice needed for the asymmetric rigid rotor, essentially to replace 
the usual Cartan generator of rotations about a fixed axis (corresponding to the magnetic quantum 
number) by a quadratic, but non-Casimir, in the generators, fortuitously turns out to belong 
technically to what in the general case is a so-called Bethe subalgebra of the associated Yang- 
Baxter algebra - in the two-dimensional case, just a function of the independent pieces of the 
transfer operator as we have seen. However, the situation becomes less fortuitous when compared 
with the analogous labelling problem for the three dimensional unitary group SU(3) when, in many 
physical applications, the natural angular momentum embedding SU(3) D SO(3) requires the use 
of the orbital angular momentum and orbital magnetic quantum numbers as usual, but also, an 
additional so-called missing label needed to remove degeneracy in the reduction of representations 
of SU (3) to their angular momentum constituents. It has been shown in [TO] that, of the known 
non-subgroup, SO(3) invariant operators in the enveloping algebra of SU{3), there is a unique 
combination which belongs to the Bethe (maximal abelian) subalgebra of the relevant Yang-Baxter 
algebra for this case. Although the spectra of the candidate operators in this SU(3) D SO(3) case 
are known to be irrational, Bethe Ansatz considerations would suggest that the eigenvalues of the 
distinguished combination are at least algebraic, and indeed, an operator derived from a Bethe 
subalgebra would appear to be best choice for such non-subgroup labelling situations - just as in 
the present problem of the asymmetric rigid rotor. 

The rigid rotor as a dynamical system is a seminal problem and has been the subject of 
innumerable papers. Its quantum version, especially the symmetric case, is a textbook example, 
while the asymmetric case is surprisingly rich, and is a basic starting point for molecular studies. 
It is of great interest that a range of new tools is available for its analysis via the techniques of 
exactly-solvable models. We defer the detailed investigation of such further implications of our 
transcription to future work. 
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